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Abstract 



First information on the chirally odd twist-3 proton distribution function e a (x) is extracted 
from the azimuthal asymmetry, Am, in the electro-production of pions from deeply inelastic 
scattering of longitudinally polarized electrons off unpolarizcd protons, which has been recently 
measured by CLAS collaboration. Furthermore parameter-free predictions are made for az- 
\ imuthal asymmetries, Aul, from scattering of an unpolarized beam on a polarized proton target 

7—1 ■ for CLAS kinematics. 

00 
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■ 1 Introduction 

O 

\ Experimental information on the chirally odd twist-3 proton distribution function e a (x) pi ^] from 

deeply inelastic scattering (DIS) would provide not only insights into the twist-3 nucleon structure. 
Qh. The first moment of e a (x) is related to the pion- nucleon cr-term, which in turn is related to the 
strangeness content of the nucleon. Here one faces the so called "cr-term puzzle". Results from 
chiral perturbation theory and the value a — (60 — 80) MeV extracted from pion-nucleon scattering 
data ||, U] imply that around 10% of the nucleon mass is due to the strange quark. This contrasts 
the fact that strange quarks carry a negligible fraction of the nucleon momentum at say 1 GeV 2 , 
the "typical hadronic scale" for nucleon set by the nucleon mass M N . 

Since e a (x) is a spin-average distribution, it can be accessed in experiments with unpolarized 
nucleons. However, due its chiral-odd nature and twist-3 character it can enter an observable only 
in connection with another chirally odd distribution or fragmentation function, and with a power 
suppression M N /Q, where Q is the hard scale of the process. So one is lead to study processes at 
moderate Q, to which e a (x) gives the leading contribution. 

An observable, where e a {x) appears as leading contribution, is the azimuthal asymmetry Am 
in pion electro-production from semi-inclusive DIS of polarized electrons off unpolarized protons 
H IHQ In this quantity e a (x) appears in connection with the chirally and T-odd twist-2 "Collins" 
fragmentation function H^- a (z), which describes the left-right asymmetry in fragmentation of a 
transversely polarized quark of flavour a into a hadron ||, ||, 0]. In the HERMES experiment 
Am was found consistent with zero within error bars M, 0], More recently, however, the CLAS 



collaboration reported the measurement of a non-zero Am in a different kinematics [10] 



So the CLAS data [10] allow - under the assumption of factorization - an extraction of first ex- 
perimental information on e a (x) from DIS, provided one knows H^. First experimental indications 



to Hi came from studies of e + e -annihilation [11]. The HERMES data on azimuthal asymmetries 



1 In Axy X[Y) denotes beam (target) polarization, and one should take the values U for unpolarized, L for 
longitudinally polarized. We use the notation of [H, H], with H^(z) normalized to (Ph±_) instead of M^. 
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Am, in pion electro-production from DIS ]8|, (9| provide further information on H^(z). In these 
asymmetries H^-(z) enters in combination with the chirally odd twist- 2 nucleon transversity dis- 
tribution hf(x) |l|, H|, [TJ|], the twist-3 distribution h a L (x) Jl], ||], and quark transverse momentum 
weighted moments thereof Q. In Ref. Q H^(z) has been extracted from the HERMES data ||, 
using for h\{x) and h%(x) predictions from the chiral quark soliton model [14] and the instanton 
model of the QCD-vacuum (if 



In this note we will use the information on H^(z) obtained in Ref. [13] to extract the twist-3 
distribution e a (x) from the CLAS data [TO]. Furthermore, we will predict azimuthal asymmetries 
Ajjl for CLAS, which are under current study. 

2 The twist-3 distribution function e a (x) 

The twist-3 quark and antiquark distribution functions e q {x) and e^(x) are defined as g 

= ^fj^ e iX *(N\M0)M^n)\N) , e«» = e«(-x) (1) 
where the insertion of the gauge-link is understood. The Q 2 -evolution has been studied in Refs. [16, 



17, [lq]. In the multi-colour limit the evolution of e a (x) simplifies to a DGLAP-type evolution - as 
it does for the other two proton twist-3 distributions h a L {x) and g^(x). The latter give a constraint 
on e a (x), the "twist-3 Soffer inequality", as follows from JTJI] 



e a {x)>2\g a T {x)\-hl{x). (2) 

At small x it behaves as, with some constants 

e a (x) x ^c l x~ 0M + c 2 d(x) . (3) 

The first term follows from Regge phenomenology e(x) ~ x~( a+1 \ However the Pomeron residue 
is, as is known, non-spin- flip, and thus decouples from the chirally odd e a (x). Therefore the small 
x-behaviour of e a (x) is determined by the lowest lying spin flip trajectory, i.e. the one with the 
scalar meson / (980). With the usual slope a' w 1 GeV" 2 this yields a rise like x ° 04 . The constant 
ci in Eq. (||) is proportional to m q jM^ due to Eq. (0) below. The second term in Eq. (|3|), the 
possibility of a <5-function at x = 0, has been recently discussed in Ref. [20]. 
The first moment of (e u + e d )(x) is related to the pion-nucleon sigma-term 

1 dx (e u + e d )(x) = -J-(iV| + ^ d ) \N) = — , (4) 

m av /AT| , T , t , \ lAn f (64 ± 8) MeV Ref. [3] 
° = 2Mn" ^ + ^ |iV) = { (79 ± 7) MeV Ref. Q . < 5 > 

With the average mass of the light quarks m av = \ (m u + m^) ~ 5 MeV one obtains 

J 1 dx (e u + e d )(x) ~ (12- 16) . (6) 

However, considering Eq. (||), this does not necessarily imply that (e u + e d )(x) itself is large. The 
second moment is proportional to the number of the respective valence quarks N q (for proton 
N u = 2 and = 1) and vanishes in the chiral limit Q 

j\xxe\x) = ^N q . (7) 
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A model estimate for quark distributions e q {x) has been given in the framework of the bag 
model ^TJ. At the estimated model scale of about 0.4 GeV the saturation of the "twist-3 Soffer 
inequality" Eq. (0) as e(x) = 2gx{x) — hz(x) has been observed pi]]. The flavour index is dropped, 



since the quark distributions of Refs. ||, 21] are flavour independent. 

Finally we remark that the twist-3 quark distribution e q (x) and the unpolarized twist- 2 quark 
distribution ff(x) coincide in the non-relativistic limit 

lim e q (x) = lim f?(x) = N q S[ x - ^ , (8) 

non non V Q / 

rclativistic rclativistic 

in which the current quark mass in Eq. (]?]) is to be interpreted as the "constituent quark" mass 
m q = | M N . The sum rule Eq. (||) is however strongly underestimated in this limit. 

3 The Collins fragmentation function 

The crucial ingredient for the extraction of the twist-3 distribution function e a (x) from the az- 
imuthal asymmetry Am is the knowledge of H^(z). 

This fragmentation function is responsible for a specific azimuthal asymmetry of a hadron in 
a jet around the axis in direction of the second hadron in the opposite jet due to transversal spin 
correlation of q and q. It was the measurement of this asymmetry, using the DELPHI data collection 



[11], which provided first experimental indication to H^. For the leading particles in each jet of two- 
jet events, averaged over z and kj_ and over quark flavours, a "most reliable" (because less sensitive 
to the unestimated systematic error) value of the analyzing power of \(H^-} / (Di}\ = (6.3 ± 2.0)% 
was found. Using the whole available range of the azimuthal angle (and thus a larger statistics) the 
"more optimistic" (and also more sensitive to the systematic errors) value for the analyzing power 



(Di) 



(12.5 ± 1.4)% [DELPHI, extraction] (9) 



was found. The result Eq. (||) refers to the scale M§ and to an average z of (z) ~ 0.4 [pj] ]. 

Combining the information Eq. (^) for with predictions for /if (x) and h°^{x) from the chiral 
quark soliton model [jl4| and the instanton model of the QCD- vacuum [|l^], it was possible to 
describe well the HERMES data on the Ajjl asymmetries & M in a parameter-free approach [ ]T~3| ] . 
For that a weak scale dependence of the analysing power Eq. (g) had to be assumed, which however 
is not supported by studies of Sudakov suppression effects [22]. 



Furthermore, in Ref. [ 13 1 — assuming the model predictions [14, 15] for the proton chiral odd 
distributions - the z-dependence of the favoured pion fragmentation function H^-(z) has been 
deduced from HERMES data []|, [|. The result refers to a scale of about 4 GeV 2 and can be 
parametrized by a simple fit 

H^(z) =azD 1 {z) with a = 0.33 ± 0.06 for 0.2 < z < 0.7 , (10) 

^— if = (13.8 ± 2.8)% for (z) = 0.41 [HERMES, extraction] , (11) 

where D\(z) is the favoured unpolarized pion fragmentation function. The errors in Eqs. (]To| ) are 
due to experimental error of HERMES data [||, ^] . The assumption of the predictions from []l4|, |l5| 
introduces a model dependence, which can be viewed as a "systematic error" and is estimated to 
be around 20%. The z-averaged value Eq. ( |i~l| ) is close to the DELPHI result Eq. (||), indicating 
that the ratio (H^) / (D\) might indeed depend on scale only weakly. Note also, that HERMES 
data favour clearly a positive sign for the analyzing power. It is noteworthy that a similar relation 
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between the favoured fragmentation functions H^~(z) and D\(z) (even close numerically!) was 
found in a recent model calculation 

In order to estimate the analyzing power for the CLAS experiment we assume the result Eq. (|l0|) 
to be valid up to z < 0.8, and to be only weakly scale dependent between HERMES (Q 2 ) = 4GeV 2 
and CLAS (Q 2 ) = 1.5 GeV 2 @. Due to the particular fit Eq. (0), the analyzing power is related 
to the average z of the experiment by {H^) / {D\) = a (z) with the constant a from Eq. (|To|). For 



the CLAS experiment [1C] we obtain 



(Hi 



(20 ± 



4 The azimuthal asymmetry A LU 



for (z) = 0.61 [CLAS, prediction]. 



sin <f> 



(12) 




i n the CLAS experiment. In the CLAS experiment a 
longitudinally polarized 4.3 GeV electron beam was scattered off an 
unpolarized proton target. The cross sections cr^ for the process 
ep — > e'7r + X were measured in dependence of the azimuthal angle 
(f), i.e. the angle between the lepton scattering plane and the plane 
defined by the momentum q of the virtual photon and momentum 
Ph of the produced pion, see Fig.l. The signs ^ refer to the 
longitudinal polarization of the electrons, with ( + ) if polarization Figure 1: Kinematics of the pro- 
parallel to beam direction, and (~) if anti-parallel. Let P (Ph) be cess ep — > e'hX in the lab frame. 
the momentum of the incoming proton (outgoing pion) and I (l 1 ) 

the momentum of the incoming (outgoing) electron. The relevant kinematical variables are center 
of mass energy square s := (P + 1) 2 , four momentum transfer q := I — I' with Q 2 := —q 2 , invariant 
mass of the photon-proton system W 2 := (P + q) 2 , and x, y and z defined by 



2Pq 



y 



2Pq 



z := 



PPh 
Pq 



(13) 



In this notation the azimuthal asymmetry A s ^j (x) measured by CLAS is given by 

/dydrd^in^ ( 1 dV+) 1 dVH ^ 

J \Si + ^ dxdydzdcj) 5^") dxdydzi 



*sin0 1 



- I dfl d : ( 



' dV+j 

dx dy dz c 



+ 



dV-> N 

dx dydzdc^) 



(14) 



where Se 1 * 1 ^ denotes the electron polarization. When integrating over y and z the experimental cuts 
have to be considered [ 10 1 



0.15 < x < 0.4 , 0.5 < y < 0.85 , 0.5 < z < 0.8 , 
1.0 < Q 2 /GeV 2 < 3.0 , 2.0 < W/GeV < 2.6 . 



(15) 



in theory. The cross sections entering the asymmetry A s ^j Eq. (14) have been computed 
in Ref.Q at tree-level up to order 1/Q. Assuming a Gaussian distribution of quark transverse 
momenta one obtains for the A 5 ^^ asymmetry Eq. (]i~4l ) 



.amy 



(x) 



JdyAy^T 



yM N /Q5j2 a e 2 a xe a (x){Ht a ) 



(z)Ji + (P 2 N }/(k 2 ) Jdy (i + (i - v) 2 ) /Q 4 n e 2 b fi(x)W, 



(16) 
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where (Ij_ 2 N ) denotes the mean square transverse momentum of quarks in the nucleon and (k]_] 



of 



the fragmenting quarks. The latter is related to the transverse momentum of the produced pion 
byfl (kl) = (P 2 hL )/(z 2 ). In the CLAS experiment (P h± ) = 0.44 GeV » (R N ) |0|. 

Eq. ( |l~4| ) assumes factorization to hold, and for that a large Q 2 is a necessary condition. Aside 
the general problem of factorization of pr-dependent processes there is a subtle question is whether 
Eq. ( |l4|) can be applied to analyze the CLAS experiment where (Q 2 ) = 1.5 GeV 2 ||i"0|| . Here we 
will assume that this can be done. This assumption will receive a certain justification, if our 
predictions on the asymmetries Am (see next section) will agree well with future CLAS data taken 
at comparably low (Q 2 ). However, one will not have a more definite answer on that, until future 
experiments performed at higher Q 2 will have constrained e a (x) such that a comparison between 
results at the different scales - taking Q 2 -evolution into account - will be possible. 



The extraction of e a (x) from preliminary CLAS data. 

Using isospin symmetry and favoured flavour fragmentation 



Z>, = D 



D 



~ 



(17) 



and the same relations for , in the expression for the az- 
imuthal asymmetry A^y Eq. (|l6|), we see that the CLAS 
data yield information on the flavour combination 



e(x) 

5 - 
4 
3 
2 



■ e(x) - 



bound \ 



1 l 



e(x) = e u (x) + ^e J (x) . 



(18) 



With the estimate of the analyzing power Eq. (12) and us- 



HERMESstat + syst 



ing for /f (x) the parameterization of Ref. [24] we obtain the 



0.1 0.2 0.3 0.4 



0.5 



Figure 2: The Savour combination e(x) = 
(e u + je d )(x) extracted from preliminary 
CLAS data vs. x at (Q 2 ) = 1.5 GeV 2 . The 
error bars are due to statistical error of 
the data. For comparison the same flavour 
combinations of /" (x) and the twist-3 Sof- 
fer lower bound are shown. 



result for e(x) of Eq. (18) shown in Fig. 2. For comparison 
the corresponding flavour combinations of the twist-3 Soffer 
bound of Eq. (||)f] and the unpolarized distribution function 
ff (x) are plotted in Fig. 2. 

Note that the uncertainties of Hi(z) in Eq. (|To|) - due 
to experimental error of HERMES data and theoretical as- 
sumptions in their analysis - affect the overall normalization 
of the extracted e(x). Its x-dependence, however, is entirely 
due to the CLAS data. 

The extracted e(x) is clearly larger than our estimate of its twist-3 Soffer bound Eq. (|2|), about 
two times smaller than ff(x) at the scale of 1.5 GeV 2 . The result indicates also that the large 
number in the sum rule Eq. (||) may require a significant contribution from the small x-region, 
which is interesting in the light of the predictions in Eq. ([|). 

It is worthwhile mentioning that the bag model result for e(x) of Ref. |2l|] (evolved according 
to naive power counting to the comparable scale of Q 2 = 1 GeV 2 ) is in qualitative agreement with 
the extracted e(x). 

2 Whether these relations hold exactly or only approximately, depends on the chosen jet selection scheme, as 
does the question, whether (kj) is a function of z. Considering the large uncertainties on both experimental and 
theoretical side, a discussion of jet selection scheme dependence seems not appropriate here. 

3 We use the "Wandzura-Wilczek approximations" g^(x) = J 1 d£ and h a L (x) = 2x £ d£h1 (£)/£ 2 - The 

neglect of the pure twist-3 h1(x) and Qt{ x ) is justified in the instanton QCD vacuum model [ fujl p5~| . For hf(x) the 
model prediction mm is used, for gf(x) the parameterization of Ref. [p6|. 
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A s ^j^ in the HERMES experiment. In the HERMES experiment the asymmetry A s ^f has 
been measured with a longitudinally polarized 27.6 GeV positron beam in the kinematical range 

0.023 < x < 0.4 , 0.2 < y < 0.85 , 0.2 < z < 0.7 , 

1 < Q 2 /GeV 2 < 15 , 2 < W/GeV , (19) 

and the following, consistent with zero result for the totally integrated asymmetries found ||] 
^V + )hermes = -0.005 ± 0.008 ± 0.004 

^^(7r-) HEHMBS = -0.007 ±0.010 ±0.004 [HERMES]. (20) 

In order to see that the CLAS [|l(| and HERMES || data are compatible we very roughly 'pa- 
rameterize' e a {x) ~ fi (x) at (Q 2 ) = 1.5 GeV 2 . This estimate is consistent with CLAS data 
(for the flavour combination (e u + j e d )(x), see Fig. 2) and describes e a (x) sufficiently well for our 
purposes. We can assume this parameterization to be valid also at the scales in the HERMES 
experiment, since evolution effects are small compared to the crudeness of our 'parameterization'. 
This allows us to estimate A s ^'(tt + ) ~ 0.008 and ^^(vr - ) ~ 0.007 for HERMES kinematics, 
which is in agreement with the data in Eq. (p0|). We conclude that the e a (x) extracted from the 
CLAS experiment (Fig. 2) is not in contradiction with HERMES data ||. 



5 Predictions for Ajjl asymmetries at CLAS 

In the HERMES experiment the azimuthal asymmetries A^f and A^^ in the production of 
charged § and neutral || pions from a proton target have been measured as functions of x and z. 
For 7r + and ir° sizeable Aff^ asymmetries have been observed, while the other asymmetries have 
been found consistent with zero within error bars. In Ref. @ the HERMES data @, |] has been 



well described in a parameter-free approach, using for the DELPHI result 0, see Eq. ®, 
and for proton transversity distributions the predictions from the chiral quark soliton model fll4|l 
and the instanton model of the QCD vacuum |]15[. This approach has been used in Ref. p7[ to 



predict Ajjl azimuthal asymmetries for a deuterium target. Here we predict A^jf and A^j for 



pion production from a proton target for CLAS in an approach similar to Ref. ]l3j, relying on 
the assumption that factorization holds at the energies of the CLAS experiment. For the CLAS 



kinematics, however, the DELPHI result [11| for H± Eq. (H) cannot be used, as it refers to a 
different (z). Instead we use our estimate from Eq. (|i~2|). Our predictions^ are shown in Fig. 3, for 
beam energies of 4.25 GeV, 5.7 GeV and 12 GeV which are currently available or proposed for the 
CLAS experiment. 

Fig. 3 demonstrates that the predicted CLAS asymmetries are as large as the asymmetries 
measured by HERMES (8|, ||]. Thus, with the high luminosity of the CLAS experiment, a precise 
measurement A s ^jf and A 3 ^ 2 ^ for ir + and 7T° is probably possible. Moreover, the CLAS kinematics 
for the 12 GeV beam allows to observe the change of sign of the A s ^ n L ^(x) asymmetries at x ~ 0.5. 
This change of sign is due to different signs of the twist-3 and twist-2 contributions. For the 
5.7 GeV beam the A s ^ n I f(x) become zero close to the upper x-cut, which makes this phenomenon 
more difficult to observe. For HERMES kinematics the zero of A^f(x) lies outside the covered 



a>range and is invisible [13|, |27fl. 



The A^jf{x) asymmetries for different pions cross each other in a single point, see Fig. 3. 
This interesting observation is due to the fact, that only two of the three cross sections for the 
production of 7r + , 7r° and tt~ are "linearly independent" because of isospin symmetry and favoured 



*For explicit expressions for the azimuthal asymmetries and further details see Refs.jl^, |2?], 
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Figure 3: Predictions for azimuthal asymmetries A v ^ (x) vs. x for different beam energies and the corresponding 
kinematical cuts at CLAS. The thick lines correspond to W(<j>) = sin</>, the thin lines correspond to W{4>) = sin 20. 
Hereby solid lines refer to 7r + , long-dashed lines to tv°, and short-dashed lines to tv~ . 



flavour fragmentation. Thus, if two curves cross each other in some point, the third one necessarily 
goes through this point as well. The exact positions of this point and of the zero of A s ^ I f{x) depend 
on the beam energy and move to smaller x with the energy growth. The experimental check of 
this prediction, especially at COMPASS energies, would give an argument in favour of the handbag 
mechanism of the asymmetry with different signs of twist- 2 and twist-3 contributions. 

Our predictions are based on the assumption that factorization holds at the scales 1 GeV 2 < 
Q 2 < 9 GeV 2 covered in CLAS experiment |1C]. It will be exciting to learn from the comparison of 
these predictions to future CLAS data, to which extent factorization holds. In particular, this will 
give valuable indications on the correct interpretation of the data on the A m asymmetry and the 
extraction of the twist-3 distribution function e a (x) given in the previous section. 

6 Conclusions 



We have presented the extraction of first information of the chirally odd proton twist-3 distribution 
function e a (x) from the azimuthal asymmetry Ajjj in 7r + electro-production from semi- inclusive 
DIS of polarized electrons off unpolarized protons, which has been recently measured by CLAS. 
The flavour combination (e u + \e d ){x) extracted in the x-region 0.15 < x < 0.4 refers to a scale 
of 1.5 GeV 2 and is sizeable - roughly half the magnitude of the unpolarized distribution function 
at that scale. But it is not large enough to explain the large number for the first moment of 
(e u + e d )(x), related to the pion nucleon sigma term, by contributions from valence x-regions alone. 

The extraction relies on the assumption of factorization, which might be questioned at the 
Q 2 of the CLAS experiment. To test this assumption, we have predicted azimuthal asymmetries 
Aul in pion electro-production from DIS of unpolarized electrons off polarized protons for CLAS 
kinematics, which are under current study. The predictions are based on a parameter-free approach, 
which has been shown to describe well the corresponding data from the HERMES experiment. A 
successful comparison of these predictions to future CLAS data would support the assumption of 
applicability of factorization at the moderate scale. 

For a definite clarification of the question, whether the CLAS data has been interpreted here 
correctly, we have to wait for data from future high luminosity (needed to resolve the twist-3 
effect) experiments performed at scales where factorization is less questioned. Maybe COMPASS 
experiment at CERN could be one of them. Our predictions for COMPASS will be published 
elsewhere. 
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